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Abstract 



Unstable states of the eta meson and the ^He nucleus predicted 
using the time delay method were found to be in agreement with a 
J> ' recent claim of r/-mesic ^He states made by the TAPS collaboration. 

^ ■ Here, we extend this method to a speculative study of the unstable 

, states occurring in the ryd and ry^He elastic scattering. The T-matrix 

I for r/^He scattering is evaluated within the Finite Rank Approxima- 

tion (FRA) of few body equations. For the evaluation of time delay in 
, the rjd case, we use a parameterization of an existing Faddeev calcu- 

lation and compare the results with those obtained from FRA. With 
an rjN scattering length, o^^r = (0.42, 0.34) fm, we find an rjd unsta- 
ble bound state around —16 MeV, within the Faddeev calculation. A 
^ I similar state within the FRA is found for a low value of o^tv; namely, 

^ I QriN = (0.28,0.19) fm. The existence of an ?7^He unstable bound state 

close to threshold is hinted by o^at = (0.28, 0.19) fm, but is ruled out 
^ , by large scattering lengths. 



PACS numbers: 14.40.Aq, 03.65.Nk,14.20.Gk 
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1 Introduction 



More than a decade after the prediction of the r^-mesic nuclei [1], their exis- 
tence remains a matter of interesting debate among nuclear physicists. The 
root of the debate lies in the uncertainty in the knowledge of the elementary 
eta-nucleon (rjN) interaction. Though once again, the first prediction of an 
attractive rjN interaction [2] (arising basically due to the proximity of the rjN 
threshold to the A^*(1535) resonance) was made back in 1985, an agreement 
on the magnitude of the attraction has not been reached. Though most of 
the models are constrained by the same sets of data on tiN elastic scattering 
and the ttN — > rjN cross sections, the rjN scattering length predicted by dif- 
ferent models is very different due to the unavailability of direct experimental 
information on rjN elastic scattering. Consequently, the predictions of possi- 
ble resonances or unstable bound states of rj mesons and nuclei within these 
models also vary a lot. In the present article we refer to 77- nucleus states 
with a negative binding energy but a finite lifetime (width) as "unstable 
bound states" . These states are sometimes also called "quasibound" states 
in literature. 

Even if the present knowledge of the rjN interaction is somewhat poor, 
the growing experimental efforts in the past few years could improve the 
understanding in a not so distant future. The photoproduction of ?7-mesic ^He 
investigated by the TAPS collaboration [3], has indeed proved to be a step 
forward in this direction. The total inclusive cross section for the 7 ^He — > rjX 
reaction was measured at the Mainz Microtron accelerator facility using the 
TAPS calorimeter and an unstable bound state with a binding energy of 
—4.4 ± 4.2 MeV was reported. Certain evidence for the existence of the 
?7-mesic nucleus ^^C, was reported in [4] and a claim for the light rj'^B.c quasi- 
bound state was made in [5] from the study of the cross section and tensor 
analysing power in the dd — > ^Hcr] reaction. Indirect evidence of the strong 
?7-nucleus attraction was also obtained from data on rj production in the 
pd, ^ ^Her] [6] and np drj [7] reactions which display large enhancements 
in the cross sections near threshold. More information is expected to be 
available from the ongoing program of the COSY-GEM collaboration [8] . 

Since the 77^He states which we recently predicted [9] from time delay were 
found to be in agreement with the TAPS data, we thought it worthwhile to 
extend the calculations and search for eta-mesic states in two other light nu- 
clei. In contrast to conventional methods of resonance extraction like Argand 
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diagrams and poles of the S'-matrix in the complex energy plane, the time de- 
lay method has a physical meaning which was noticed more than 50 years ago 
by Eisenbud [10] and Wigner [11]. The fact that a large positive time delay 
in a scattering process is associated with the formation of a resonance is now 
mentioned and elaborately discussed in most standard text books [12]. This 
time delay is related to the energy derivative of scattering phase shifts and 
can be easily calculated. However, strangely but not surprisingly, this well- 
documented method had rarely been used in literature to extract information 
on resonances from data until recently when the authors of the present work 
put it to a test with hadron-hadron elastic scattering [13, 14, 15, 16]. Be- 
ing encouraged by the fact that this method was successful in characterizing 
known nucleonic resonances [16], could find some evidence for pent a- quark 
baryons [13] and revealed all known meson resonances [14] , we carried out a 
similar study [9] for the ry^Hc system which also proved quite fruitful. Hence, 
in the present work we shall search for r]d and T^^He unstable states using the 
time delay method. The transition matrix for ^^-nucleus scattering (which is 
required as an input to the calculation of time delay) is constructed using few 
body equations and accounts for off-shell re-scattering and nuclear binding 
energy effects. This model was successful [17] in showing that the enhance- 
ment in the cross sections (or the strong energy dependence of the scattering 
amplitude near threshold) of the pd — > ^He rj reaction is due to the rj '^He 
final state interaction. Since the calculations are purely theoretical, we have 
extended the concept of time delay to negative energies too. In the following 
sections, we shall demonstrate the validity and usefulness of this concept in 
characterizing bound and virtual states occurring at negative energies. 

There exist scores of papers in literature, with predictions of unstable 
bound (sometimes called quasi-bound) states of ?7-mesic nuclei as well as rj- 
nucleus resonances. One can find a detailed account of the existing literature 
in a recent work [18], where a search for unstable r^-mesic nuclei ranging from 
A = 3 to 208 was made. Here, we shall briefly survey the literature on light 
nuclei, since we restrict ourselves to the study of ^^d and r^^He elastic scatter- 
ing in the present work. One of the early calculations [19] using the 3-body 
equation, predicted a resonance with a mass of 2430 MeV and a width of 10- 
20 MeV in the TrNN-rjNN coupled system. This state led to a remarkable 
enhancement of the rjd elastic cross section. In some other Faddeev-type cal- 
culations [20, 21] of r]d scattering, a resonance at low energies was predicted. 
The exact values of the predicted mass and the width of course varied with 
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the value of the r^-nucleon scattering length, a^TVj which in turn depends on 
the model of the ?7-nucleon interaction. The existence of resonances in the 
above works was inferred from Argand diagrams. Some other recent (also 
Faddeev-type) calculations, however, ruled out the possibility of a resonance 
in the rjd system [22, 23, 24]. Predictions of resonances in the light ?7-nucleus 
systems, namely, d, "^Hc and ''He, from the positions and movements of the 
poles of the amplitude, can be found in [25]. Finally to mention some pre- 
dictions of 'virtual' states, a virtual (anti-bound) s-wave rj^R state, which 
led to a large enhancement of the cross section for r] production from the 
three-body nuclei was found in [26], whereas a narrow virtual state in the 
Tjd system was found [27] to have a rather weak effect on the pn — > rjd cross 
sections. 

In this work, we shall infer the existence of unstable states in the rjd and 
r/^Hc systems, from positive peaks in time delay and compare our results with 
existing predictions in literature using other methods. In the next section, we 
present the basic elements of the "time delay" method and the characteristics 
expected from it in the case of resonances, bound, quasi-bound, virtual and 
quasi-virtual states. Since one can find thorough discussions of time delay 
in literature [28, 29], and also in standard text-books on quantum mechanics 
[12, 30], we do not perform a review of this method here. However, the 
efficacy of the inferences from this method is established by applying it to 
the known case of neutron-proton bound and virtual states and referring to 
our earlier application to the known hadron resonances. Section 3 is devoted 
to a brief discussion of the few body equations involving the finite rank 
approximation (FRA) which we use to calculate the ?7-nucleus transition 
matrix. This t-matrix is subsequently used to evaluate the time delay in 
?7-light nucleus scattering. Since the applicability of the FRA for the rjd 
case is fimited (it was shown in [21] that FRA agrees with the Faddeev 
calculations for real part of the scattering lengths less than 0.5 fm only), 
we use a parameterization of the t] d t-matrix using relativistic Faddeev 
equations and present the results of this calculation in Sec. 4. The time 
delay results for the T^^'He system are given in Sec. 5. The FRA calculations 
are done using two different models of the rjN interaction (constructed within 
coupled channel formalisms) which fit the same set of data on pion scattering 
and pion induced t] production on a nucleon, but give different values of the 
rjN scattering length. Finally, Sec. 6 summarizes the findings of the present 
work. 
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2 Time delay plots of bound, virtual and de- 
caying states 

The collision time or delay time in scattering processes, was quantified by 
Eisenbud and Wigner in terms of the observable phase shift which can be 
extracted from cross section data. In [11], Wigner pictured the resonance 
formation in elastic scattering as the capture and retention of the incident 
particle for some time by the scattering centre, introducing thereby a time 
delay in the emergence of the outgoing particles. He further showed that the 
energy derivative of the phase shift, S, which is related to the time delay, 
At{E), as 

At{E) = 2h^ (1) 

and is large and positive close to resonances, can also take negative values 
which are however limited from causality constraints. In the presence of 
inelasticities, a one to one correspondence between time delay and the lifetime 
of a resonance does not hold and a more useful definition, namely, the time 
delay matrix (later discussed in terms of a lifetime matrix by Smith [29]) was 
given by Eisenbud. An element of this matrix, Atij, which is the time delay 
in the emergence of a particle in the j^^ channel after being injected in the 
i*'* channel is given by. 



Atij = 



(2) 



where Sij is an element of the corresponding ^'-matrix. Writing the S'-matrix 
in terms of the T-matrix as, 

S = l + 2iT, (3) 

one can evaluate time delay in terms of the T-matrix. The time delay in 
elastic scattering, i.e. Ata, is given in terms of the T-matrix as. 



(4) 



SlS,,At^, - 2% , „ --ydEj " -^KdE 

where T is the complex T-matrix such that, 

Tfcj = ^eTkj + i ^mTkj . (5) 
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The time delay evaluated from eqs (1) and (4) is just the same. 

The above relations were put to a test in [13, 14, 15, 16] to characterize 
the hadron resonances occurring in meson-nucleon and meson-meson elastic 
scattering. The energy distribution of the time delay evaluated in these 
works, nicely displayed the known and A baryons, meson resonances like 
the p, the scalars (/o) and strange i^'*'s found in Kti scattering, in addition 
to confirming some old claims of exotic states. The peaks in time delay, 
At{E), agreed well with the known resonance masses. It was also shown 
that the time delay peaks and the T-matrix poles essentially contain the 
same information. A theoretical discussion on this issue can be found in 
[30, 31]. 

The time delay concept is not only useful to locate resonances, but can 
also be used to locate the bound, virtual and unstable bound states which 
have negative binding energies. We illustrate this assertion with the well- 
known case of the n—p system. The S'-matrix for the neutron-proton system, 
constructed from a square well potential which produces the correct binding 
energy of the deuteron is given as a function of / as, 

g ^ ahf{a)m-(3hf\a)j{{(3) 
' ahf{a)m-Ph?{a)m 

where j;, hi'p and hf^ are the spherical Bessel and Hankel functions of the 
first and second kind respectively. 

a^kR and p = {a^ - 2iiU /h'^f'^ (7) 

where the potential U is given by 

U = V + iW, U{r) = Ue{R - r) (8) 

and R is the width of the potential well. A similar square well potential 
was used by Morimatsu and Yazaki [32] while locating the "unstable bound 
states" of S-hypernuclei as second quadrant poles and by J. Fraxedas and J. 
Sesma using the time delay method [33]. We evaluate the time delay in n — p 
scattering at negative energies E, where E — \/s—mn—mp with ^Js being the 
energy available in the np centre of mass system. Using the above S'-matrix 
with a = ikR (hence E = —k'^/2fi), I = and the appropriate parameters 
for an n — p square well potential, namely, V = Si.QMeV and R = 2.07 fm. 
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Figure 1: The theoretical delay tmie in np clastic scattering as a function of the 
energy E = ^/s — m,n — rup, where -s/s is the total energy available in the np centre 
of mass system. 



the time delay plot (as in Fig. 1) shows a sharp spike (positive infinite time 
delay) exactly at the binding energy of the deuteron {E = —2.224 MeV). 
If we add a small imaginary part to the potential, then of course there is 
a Breit-Wigner kind of distribution centered around the binding energy of 
the deuteron (a fictitious "unstable bound state" of the n — p system at 
E = -2.224 MeV). 

The correlation between the potential parameters and the position of the 
spike at the correct deuteron binding energy is very definite. If we take 
the potential parameters which do not give the correct binding energy, then 
the spike appears at a wrong place in the time delay plot. In Fig. 2, we 
demonstrate this sensitivity of time delay. In the upper half of the figure, 
we plot time delay calculated using a fixed well depth of = 34.535 MeV 
and different choices for the width of the square well. It can be seen that 
the position of the spike is very sensitive to the value of R. The spike at 
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the right binding energy is produced only with R = 2.07 fm. In the lower 
half of the plot, we perform a similar calculation but this time with the well 
width fixed and the well depth changing. Again, a small change in the well 
depth parameter, changes the energy at which the positive infinite time delay 
appears. 
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Figure 2: Sensitivity of time delay in np scattering to the square well parameters. 
V is the depth of the square well potential and R its width. The solid lines indicate 
the sharp positive infinite time delay at the correct deuteron binding energy. 

Beyond quasibound states, an S'-matrix pole in the third quadrant of 
the complex momentum plane corresponds to a quasivirtual state, which 
translates to a pole on the unphysical sheet of the complex energy plane of 
the type —\E\ +i|r/2|. This, in contrast to a resonance pole of — i|r/2| 
(which leads to an exponentially decaying state with a decay law of e~'"*), 
gives rise to an exponential growth, namely e~^^^. One can then see that in 
contrast to the time "delay" that one observes for a resonance, one would 
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observe a time "advancement" for a quasivirtual state. In other words, for 
a quasivirtual state one observes a finite 'negative' time delay. Similarly, a 
virtual state, in contrast to a bound state would show an infinite negative 
time delay. This is indeed seen in Fig. 1 for the known n — p virtual state 
at 100 keV, where the time delay calculated from the square well potential 
with parameters corresponding to this virtual state (namely, V = 23.6 MeV 
and R = 2 fm) is plotted. 

3 T-matrix for ?7-nucleus elastic scattering 

We evaluate the transition matrix for ry-nucleus (rjA) elastic scattering, using 
few body equations for the r]{2N) and r){4:N) systems. The calculation is done 
within a Finite Rank Approximation (FRA) approach, which means that in 
the intermediate state, the nucleus in rjA elastic scattering remains in its 
ground state. Since the 77-mesic bound states and resonances are basically 
low energy phenomena, it seems justified to use the FRA for calculations 
of the present work. In [21], the authors mention that though the use of 
FRA for rj^He and T^'^He systems seems justified, it is questionable for the 
case of rjd scattering and investigate the shortcoming due to the neglect of 
excitations of the nuclear ground state in 77-deuteron calculations. Within 
their model, they find that the FRA results differ from those evaluated using 
the rigorous Alt-Grassberger-Sandhas (AGS) equations in the case of the 
strong 7]N interaction (Re a^jN > 0.5 fm), while for small values of Re a^Tv, 
the FRA is reasonably good. Therefore, for the ?7-d case, we also include 
calculations using the results from the recent relativistic Faddeev equation 
(RFE) calculations for the rjNN system. 

The target Hamiltonian Ha, in the FRA is written as [36], 



where ipo is the nuclear ground state wave function and £ the binding energy. 
The TjA T-matrix in the FRA is given as [25, 36, 37], 



Ha ~ e\iJo >< ipo 



(9) 



< k';^Po\t°{z)\k; V'o > + 

f dk" < k'-il:o\ t%z) \ k" ; > 



(10) 




triA{k" , k ; z^ 
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where z = E — \e\ + iO. E is the energy associated with rjA relative motion, 
6 is the binding energy of the nucleus and fi is the reduced mass of the rjA 
system. Though the operator describes the scattering of the rj meson from 
nucleons fixed in their space position within the nucleus, it differs from the 
usual fixed center t-matrices. Here, is taken off the energy shell and in- 
volves the motion of the rj meson with respect to the center of mass of the 
target. The present scheme should not be confused with a conventional opti- 
cal potential approach which involves the impulse approximation and omits 
the re-scattering of the 77 meson from the nucleons. The matrix elements for 
i° are given as, 

< k'; 'ilJo\t^{z)\k; ^po >= j df\^Q{r)\'^ f {9 , k]f] z) (11) 

where, 

(k', k-T-z)^Y. ti (k', k-n-z) (12) 

1=1 

t° is the t-matrix for the scattering of the ?7-meson from the i*^ nucleon in 
the nucleus, with the re-scattering from the other (A-1) nucleons included. 
It is given as, 

+^ (hi h " \ — +'^^ (u h " f ^^'^ ^" ' ' ^) \^ u ^ ^ 

Ij^ [K , K ,ri , Z ) — Ij^ [K , K ,ri , Z)-\- j^772 2^ Ij [K , K ,rj , Z) 

(13) 

The t-matrix for elementary ry-nucleon scattering, t^^, is written in terms of 
the two body rjN matrix tnN^riN as, 

f^^ik', k;fl;z) ^ ^N^r^Nik'-, k ; z) exp[i{ k-k')- rl ] (14) 

The ^He nuclear wave function, required in the calculation of the T-matrix 

is taken to be of the Gaussian form. The deuteron wave function is written 
using a parametrization of the wave function [34] obtained using the Paris 
potential. The results using the Paris potential are also compared with a 
calculation using a Gaussian form of the deuteron wave function. 

As mentioned in the introduction, there exists a lot of uncertainty in 
the knowledge of the 77- nucleon interaction and hence, we use different pre- 
scriptions of the r/-N t-matrix, tj^ n~*tiN, leading to different values of the 
r]N scattering length. We give a brief description of two of these models of 
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tnN^r)N which we use for the FRA calculations below. In [24] a coupled 
channel t-matrix including the ttN and rjN channels with the Sn - rjN inter- 
action playing a dominant role was constructed. The t-matrix thus consisted 
of the meson - N* vertices and the N* propagator as given below: 

tr^N^nN{k\ k]z) = j^f^-^T^^{z)-^^^-^ (15) 

with, 

r^. {z) = {z-Mo- T.^{z) - E^(;^) + le (16) 

where T,a{z) {a = tt, rj) are the self energy contributions from the ttN and rjN 
loops. We choose the parameter set with g^, = 2.13, P — 13fm~^, and Mq — 
1656 MeV which leads to a^^v = (0.88, 0.41) fm. 

We also present results using one of the earhest calculations of the 77- 
N t-matrix [2] which gives a much smaller value of the scattering length, 
namely, a^^r = (0.28, 0.19) fm. In this model, the vrN, t^N and vrA (ttttA^) 
channels were treated in a coupled channel formalism (so that an additional 
self-energy term appears in the propagator in Eq. (16)). The parameters of 
this model are, g^, = 0.616, P = 2.36 fm~\ and Mq = 1608.1 MeV. 

There also exists a recent model of the r]N interaction [40] , which predicts 
a scattering length of flj^iv = (0.91,0.27) fm, from a fit to the nN ttN, 
ttN r]N, 'jN — ^ ttN and rjN data. However, we have not used 

it for our present FRA calculations since the T-matrix which fits the data 
very well is an on-shell T-matrix. The off-shell separable form given by the 
authors [40] agrees with their on-shell T-matrix (which fits data) but does 
not include the intermediate off shell tt and rj loops. The off shell nature 
appears only in the vertex form factors. 

The T-matrix for rjA elastic scattering, t^^, is related to the 5*- matrix as, 

5 = (17) 

TT 

where k is the momentum in the rjA centre of mass system and hence, the 
dimensionless T-matrix required in the evaluation of time delay as given in 
eq. (4) is evaluated using the relation, 

2 TT 

We shall present the time delay plots for T^d and r] "^He elastic scattering in 
the next sections. 
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4 The rj deuteron system 



We make an analytic continuation of the T matrix for r^-nucleus elastic scat- 
tering on to the complex energy plane. Evaluating the matrix elements of the 
?7-nucleus T-matrix at negative energies (corresponding to purely imaginary 
momentum), i.e. ^Aiik, ik]z), we evaluate the time delay in ?7-nucleus elas- 
tic scattering and search for the "unstable bound states" . The resonances at 
positive energies are of course determined from the positive time delay peaks 
at positive energies and real momenta. 
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Figure 3: The delay time in r/-deuteron elastic scattering as a function of the 
energy E = y/s — nirj — nid, where y/s is the total energy available in the rjd centre 
of mass system. The shaded curves are calculations using the Gaussian form of the 
deuteron wave function and dashed lines are evaluated using the Paris deuteron 
wave function. The vertical axis scale is broken in order to display the structure 
in the time delay plot clearly. 

In Fig. 3, we plot the time delay in 77 d — > d elastic scattering with 
two different inputs for the elementary rjN interaction. The rjN scattering 
lengths of ttrfN = (0.88,0.41) fm and ttrfN = (0.28,0.19) fm correspond to 
rjd scattering lengths of a^^ = (1.52, 2.57) fm and a^jd = (0.67, 0.42) fm 
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respectively. In both cases, we see a large positive time delay located near 
threshold. The solid lines (with shaded regions) are the calculations using a 
Gaussian form and the dashed lines with a Paris potential parametrization 
of the deuteron wave function. The vertical axis scale in Fig. 3 is broken in 
order to display the structure in the time delay plot clearly. We do not find 
a big difference in the results with the change of the wave function. In the 
case of the weaker rjN interaction, i.e. a^^v = (0.28,0.19) fm, there appears 
a very broad bump around —15 MeV which could be due to an unstable 
bound state. On the positive energy side, there is a sharp negative time 
delay when the kinetic energy of the rj-d system equals the binding energy 
of the deuteron. This behaviour is expected because of the connection of 
the energy derivative of the phase shift and hence the time delay (1) to 
the density of states as given by the Beth-Uhlenbeck formula [41]. It was 
shown in [15] that a maximum negative time delay occurs at the opening 
of an inelastic threshold. In the present case, the negative dip around 2.22 
MeV, corresponds to the break up threshold of the deuteron. We also see a 
resonance just near this inelastic threshold. However, the above result near 
the inelastic threshold should be taken with some caution since we have used 
the FRA which might not be a very good approximation at energies where 
new thresholds open up. 

In order to check the validity of the above FRA approach for the rjd 
case (where it is known to have limitations [21]), we evaluate time delay 
using a model [42] which obtains the r]d elastic scattering amplitude using 
a relativistic version of the Faddeev equations described in [43]. The r]d 
amplitude in [42] is parametrized using the effective range formula: 

= ^ + iRvdh"" + Sr^dk^ -ik (19) 

where k is the momentum in the r^d centre of mass system and the parameters 
\d-, Sr^d and Rrfd are as given in Table II of [42]. In Figs 4 and 5, we show 
the results obtained using the above amplitude. The time delay in Fig. 4 
has been evaluated using k —ik and hence the negative dips in this figure 
correspond to the quasivirtual states which appear as poles in the third 
quadrant of the complex k-plane. The locations of these dips (using different 
models as given in [42] ) are exactly at the energy pole values given in Table 
III of [42] as expected. As explained already in Sec. II, such a negative time 
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Figure 4: Quasivirtual states of the rjd system evaluated from a Faddeev calcu- 
lation for T]d scattering. The various model numbers correspond to the different 
strengths of the r)N interaction as explained in [42]. 
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Figure 5: Quasivirtual and quasibound states of the rjd system evaluated from a 
Faddeev calculation for rjd scattering using model as mentioned in [42]. 
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Figure 6: Magnitude of the complex amplitude 5* = 1 + 2ikf in the complex 
momentum plane evaluated from a Faddeev calculation for rjd scattering using 
model as mentioned in [42]. The sharp pole corresponds to the quasibound state 
at —16 MeV as seen in the time delay plot in Fig. 5. 

delay (or time advancement) is expected for quasivirtual states which give 
an exponential rise rather than an exponential decay law. 

In the upper half of Fig. 5, we plot the time delay evaluated using k — >■ 
—ik and the lower half shows the time delay plot corresponding to k ^ +ik 
in the amphtude with Model 0. Thus the negative dip in the upper half is the 
quasivirtual state as also given in Table II of [42] and the lower half shows 
the positive time delay corresponding to a quasibound or what we address 
as an "unstable bound" state in the present work. In order to demonstrate 
once again the one-to-one correspondence between the time delay peak and 
(S'-matrix poles, in Fig. 6 we plot the magnitude of 5" as a function of the real 
and imaginary parts of momentum k. The pole occurs at (—0.283 -|- iO.674) 
fm in the complex momentum plane. This corresponds to a peak value of 
about —17 McV and a width of 35 MeV which is in good agreement with the 
time delay plot in Fig. 5. Such a pole has however not been mentioned by 
the author in [42] (probably due to the fact that the author in [42] has been 
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mostly concerned about the effect of the near threshold quasivirtual states 
on the np — > dr] reaction). We do not find any such positive pealcs in any of 
the models shown in Fig. 4. It is interesting to note that as expected from 
[21], the FRA has some agreement with the Faddeev calculation for small 
r]N scattering lengths, o^at, and the discrepancy increases for large a^^v- 

Finally, we wish to caution the reader regarding the interpretation of time 
delay peaks in the case of s-wave scattering. To see this, substituting the 
phase shift expression, S — exp{2iS) and comparing it with (17), one can 
write, 

S = ^ ln(l - ^ V) = ^ ln(l + 2ikf) (20) 

where / is the scattering amplitude. For small k,5 kf and the behaviour of 
dS/dE (the real part of which is essentially the time delay) is determined by 
the simple pole at = (or E^fA — Ethreshoia) and the energy dependence of 
the scattering amplitude /. In the absence of a resonance, as A; 0, 5 = fca, 
where a = aR + iai is the complex scattering length. For positive energies, 
3?e5 = kttR, whereas for energies below zero, k ^ ik and 3?e5 = —kai. In such 
a situation, ^e{dS / dE) exhibits a sharp peak at Ejja — Efhreshoid, the sign of 
which is determined by the sign of the scattering length. On the other hand, if 
the scattering amplitude has a resonant behaviour near threshold, one would 
see a superposition of the two behaviours. Consequently, a state reasonably 
close to threshold gets distorted in shape and one very close manifests simply 
by broadening the threshold singularity. A state far from threshold, however, 
remains completely unaffected. 

5 The r/^He system 

The time delay plot in the upper half of Fig. 7, shows once again the 
near threshold peak and one more broad one centered around —2 MeV. 
However, the plot with time delay evaluated using the model which gives 
«??Ar = (0.88, 0.41) fm, shows a large negative time delay near threshold. The 
negative time delay can sometimes also arise due to a repulsive interaction 
[29]. Intuitively, an attractive interaction is something that causes resonance 
formation and "delays" the scattering process. A repulsive interaction on the 
other hand would speed up the process and the time taken for the process in 
the absence of interaction would be larger than that with interaction. 
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Figure 7: The delay time in ry '^He elastic scattering as a function of the energy 
E = — rrirj — m^4ffg, where s/s is the total energy available in the rj^He centre 
of mass system. 



Before ending this section, we note that the ri'^B.e scattering lengths 
obtained within the FRA and using two different models of the riN in- 
teraction, namely, Gj^n = (0.88,0.41) fm and o^at = (0.28,0.19) fm are 
O'rj'iHe — (—3.94, 5.575) fm and a^4^g = (1.678, 1.524) fm respectively. 



6 Summary 

In conclusion, we summarize the present work as: 

(i) we have made a search for the unstable states of ry-mesic deuteron and 
■^He extending the approach of time delay which was used recently for the 
first time in the eta-mesic case [9] . 

(ii) The established time delay method for searching resonances has been ex- 
tended to negative energies, to search for bound, virtual and unstable bound 
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states. The validity of this method has been estabhshed by applying it first 
to the known case of the np system and then to the case of the rj-d system 
within a parameterized Faddeev calculation. 

(iii) The calculations were performed with different values of the rjN scatter- 
ing length considered as acceptable in literature. Within the Faddeev equa- 
tion parameterization, we find one unstable bound state far from threshold 

( 16 MeV) for an rjN scattering length of (0.42, 0.34) fm. Within the FRA 

calculation, we find such an rid state around —12 MeV for a^jv = (0.28, 0.19) 
fm. These results seem to indicate that though the FRA in general is not 
recommendable for rjd elastic scattering, the results are close to those from 
the Faddeev calculations for low values of the rjN scattering length. 

(iv) In the ?7*^He case, within the FRA calculations, we find an unstable bound 
state close to threshold for a small scattering length of (0.28, 0.19) fm. 
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